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Abstract 

We analyze the phase diagram of QCD with two massless quark flavors in the space of 
temperature, T, and chemical potential of the baryon charge, /x, using available experimen- 
tal knowledge of QCD, insights gained from various models, as well as general and model 
independent arguments including continuity, universality, and thermodynamic relations. A 
random matrix model is used to describe the chiral symmetry restoration phase transition 
at finite T and /x. In agreement with general arguments, this model predicts a tricritical 
point in the T/x plane. Certain critical properties at such a point are universal and can be 
relevant to heavy ion collision experiments. 
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1 Introduction 



Current and projected experimental progress in the physics of heavy ion collisions increas- 
ingly demands better theoretical understanding of the underlying phenomena. In particular, 
the most exciting possibility offered by such experiments is the creation of high temperature 
and density conditions under which the dynamics of QCD can bring matter into a new state. 
The challenge is then to calculate the properties of this new phase together with the prop- 
erties of the phase transition from QCD, the underlying theory of quark-gluon interactions. 

Substantial progress has been achieved in our understanding of QCD at high temperature, 
T. The foundation of this understanding is provided by lattice field theory Monte Carlo 
calculations. In particular, we know that in QCD with two massless flavors a transition 
restoring chiral symmetry occurs at a temperature of approximately 160 MeV JlJ . 

On the other hand, little is known about the behavior of QCD for finite baryon charge 
density, or chemical potential of the baryon charge, /x. Standard lattice Monte Carlo tech- 
niques cannot be applied since the determinant of the Dirac operator is complex, and hence 
the Euclidean path integral defining the theory does not have a Gibbs (i.e., real, positive- 
definite) measure. A Gibbs measure is needed for the probabilistic interpretation which 
forms the basis for importance sampling methods such as Monte Carlo calculations. More- 
over, the approximation of quenched fermions fails in this case || |J for reasons which have 
been understood recently using the random matrix theory [|J. However, the conditions 
created in heavy ion collision experiments require an understanding of the regime of high 
baryon density as well as that of high temperature. As reviewed, e.g., in ||, there is good 
evidence that central part of the collisions can be described approximately before freeze-out 
by thermodynamics, so that the temperature and chemical potential can be defined. 

The purpose of this paper is to assemble available knowledge about QCD and apply it 
to the construction of the phase diagram in the T\x plane. Most of the studies of this phase 
diagram have concentrated on modeling the properties of the chiral phase transition (see, 
e.g., H, [7|, H, H, [lOfl). In this paper, we present a more complete and less model dependent 
analysis of the phase diagram which also includes effects from other phase transitions, such as 
the nuclear matter liquid-gas transition. Naturally, many of the phenomena to be discussed 
have been studied extensively. As a result, we will repeat some familiar experimental facts 
and theoretical arguments (with references to some of the original papers or reviews as 
appropriate). The aim of our analysis is to transform this knowledge into the determination 
of a phase diagram for QCD in the Tfj, plane. Such an analysis is especially important as 
an extension of Monte Carlo studies, given the technical problems that these encounter with 
finite baryon charge density. 

The chiral phase transition is of primary interest in ultrarelativistic heavy ion experiments 
since this is the transition that separates the hadronic phase from the quark-gluon phase. In 
section [| we introduce a random matrix model of the chiral phase transition at finite T and 
/i. We find that this model predicts a tricritical point in the T/j, plane in agreement with 
more generic arguments. We analyze the properties of some thermodynamic observables in 
the vicinity of this point. 
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2 Definitions 



We take as our model the standard approximation in which we (i) consider pure SU(3) 
QCD with electroweak interactions turned off and (ii) consider this theory with two massless 
quarks. There is then an exact SU(2) £ xSU(2) fi xU(l)B global symmetry of the action, which 
is spontaneously broken down to SU(2)y xU(l)s at zero and sufficiently low temperatures by 
the formation of a condensate, (ipip)- Many features of QCD indicate that this is a reasonable 
approximation, e.g., the lightness of pions, the success of current algebra relations, etc. (We 
will comment below on the inclusion of electromagnetic interactions and strange quarks.) 
This theory is described by a grand canonical partition function which, when written as a 
path integral, is formally: 



Z = e - n{T '^ /T = J VAVipVipexp{-S E } . 



The Euclidean action, Se, is given by 



rl/T r 

Se = J dxo J d 3 x 



i - 
— TtF^f^ - 2^ Vv 



(2) 



where Nf = 2 is the number of flavors, N c = 3 is the number of colors, and rrif = m = is 
the quark mass. The Euclidean matrices 7 M are hermitean. Note that with our sign choices 
positive m and \i induce positive (ipip) and (ip'-foip) . The normalization of fi differs from the 
normalization customary in lattice calculations by a factor 1/N C (i.e., the baryon charge of 
a quark). Integrating over the fermion fields we can also write: 



Z = J VA exp 



1 



2g 



TrF F 

2 A 1 ^ A*f 



det 



+ m f + A 7o 



(3) 



As indicated, this system is characterized by equilibrium values of T and \i. This may 
be thought of by imagining the system to be in thermodynamic equilibrium with a large 
reservoir of entropy and baryon charge which is characterized by these values of T and 
/i. The total energy and baryon charge of our system fluctuate. Of course, the relative 
magnitude of these fluctuations is negligible for an open system of macroscopic size. The 
relation between the chemical potential, /i, and the average baryon number density (per unit 
volume), n, is the same as that between the temperature, T, and the average entropy density 
(per unit volume), s: 



sV = - — 



on 
df ' 



(4) 



where Q is the thermodynamic potential defined in eqn. ([I]). It can also be seen that Q = 
—pV, where p is the pressure. In other words, pressure, temperature and chemical potential 
are not independent variables for our system. Their variations are related by 



dp = sdT + ndfi. 



(5) 
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Both T and p (as well as p) are intensive parameters. For a system in thermodynamic 
equilibrium, these quantities are the same for any of its smaller subsystems. In contrast, the 
extensive densities s and n can differ for two subsystems even when they are in equilibrium 
with each other. This happens in the phase coexistence region, e.g., a glass containing water 
and ice. It is more convenient to describe the phase diagram in the space of intensive pa- 
rameters, T and p. In particular, the first-order phase transition which we shall encounter 
is characterized by one value of p but two values of n — the densities of the two coexisting 
phases. Another reason for working in these coordinates is that first-principle lattice calcu- 
lations are performed in such a way that T and p are the parameters that can be controlled 
while the densities are measured. The results of relativistic heavy ion collision experiments 
are also often analyzed using this set of parameters 0. 

3 Zero Temperature 

We begin by considering the phase diagram as p is varied along the line T = 0. Strictly 
speaking, we are not dealing with thermodynamics here since the system is in its ground 
state. This fact leads to a simple property of the function n(p). Let us rewrite the partition 
function, eqn. ([!]), as the Gibbs sum over all quantum states, a, of the system: 



where each state is characterized by its energy, E a , and its baryon charge, N a . In the limit 
T — > 0, the state with the lowest value of E a — pN a makes an exponentially dominant 
contribution to the partition function. When p — 0, this is the state with N = and E = 0, 
i.e., the vacuum or a = 0. Let us introduce 



As long as /i < po, the state with the lowest of E a — pN a remains the vacuum, a = 0. 
Therefore, we conclude that at zero temperature 



What is the value of /io? As an exercise, we first consider a free theory of massive fermions 
carrying one unit of baryon charge. The states which minimize E a /N a are states with one or 
two (more if fermions have flavor or other degeneracy) fermions at rest with p = 0. For each 
of these states, E a /N a = m, the mass of the fermion. Therefore, po = m for such a theory. 
When p > m, the ground state is the Fermi sphere with radius pp = y/ p? — m 2 . Therefore, 
n(/i) = (/i 2 — m 2 ) 3 / 2 /(37r 2 ). Thus, we see that, even in a trivial theory, the function n(p) 
has a singularity at p, — /io- The existence of some singularity at the point p, — po, T — 
is a robust and model independent prediction. This follows from the fact that a singularity 
must separate two phases distinguished by an order parameter, e.g., n. The function n = 
cannot be continued to n ^ without a singularity. 

What is /io for the case of QCD, and what is the form of the singularity? The answers to 
these questions are somewhat different in QCD and in the real world (QCD+) which includes 




(6) 




(7) 



n(p) = for p < p . 
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other interactions, most notably electromagnetic interactions. Since QCD is the focus of the 
present paper and QCD+ is the ultimate goal of our understanding, we shall consider both 
cases. It is important to understand their differences if we are to extract physically useful 
predictions from lattice calculations, which are performed for QCD rather than QCD+. 

The energy per baryon, E/N, can also be written as tun — (Nm^ — E)/N, where rriN = 
m p ~ m n is the nucleon mass. Therefore, the state which minimizes E/N is that for which 
the binding energy per nucleon, e = (Nm^ — E)/N, is a maximum. Empirically, we know 
that this state is a single iron nucleus at rest with N = A = 56 and e ~ 8 MeV. However, 
in QCD without electromagnetism the binding energy per nucleon increases with A. This 
is the consequence of the saturation of nuclear forces and can be seen from the Weizsacker 
formula. Without electromagnetism, only the bulk and surface energy terms are significant 
for large A: 

4A) S Am «- m * ~ ai -a^-* (9) 



with ai ~ 16 MeV, a 2 ~ 18 MeV [11 1. As A — > oo, e saturates at the value a\. This 



corresponds to the binding energy per nucleon in a macroscopically large sample of nuclear 



matter as defined by Fetter and Walecka in ||11|| . We conclude that in QCD the density 
jumps dX n = ~ m N — 16 MeV to the value of the nuclear matter density n « 0.16 fm~ 3 . 
Therefore, in QCD there is a first-order phase transition, characterized by a discontinuity in 
the function n(//) at \i = /io (see Fig. [l|a). 

In QCD+, the Coulomb forces change the situation near fiQ. The contribution of the 
Coulomb repulsion to e(A) is negative: — (0.7 MeV) Z 2 /A 2 ^ 3 , and it is responsible for the 
experimentally observed maximum in e(A) at A ~ 56. Isospin singlet nuclear matter (A = 
oo) is unstable at zero pressure due to Coulomb repulsion. Neutron matter with Z < A 
is also unstable at zero pressure, and we are left to consider a gas of iron nuclei. In order 
to ensure electric neutrality, we must add electrons. Such a gas is clearly unstable at small 
densities and forms a solid — iron. Therefore, there is a discontinuity in the value of n(fi) at 
fiQ ~ rriN — 8 MeV. This discontinuity is equal to the density of normal matter (i.e., iron) and 
is about 10~ 14 times smaller than in QCD. For very small fi — //q, n(/i) has structure, fine on 
the scale of QCD, which reflects the properties of normal matter under pressure. Then, for 
H — fi = O(10 — 200 MeV), we traverse the domain of nuclear physics with the possibility 
for various phase transitions. In particular, a transition to neutron matter [Z <C A) is 
probably similar to the transition in QCD at // = /io- (See Fig. [l|b.) In this domain, one 
may encounter such phenomena as nuclear matter crystallization [|T^, |T3|| , superconducting 



phases of neutron and quark matter [13, O, TB\, and, due to the strange quark in QCD 



kaon condensation |12|, [17]] and a transition to strange quark matter Jl8|, Moving along 
the [A axis to the right is equivalent to increasing the pressure: p = J ndfi. Thus, this picture 
is roughly what one might encounter in moving towards the center of a neutron star from 
the iron crust at the surface. 

Our knowledge of n(/i) is scanty for densities of order one to ten times n and fi — 
fi = O(10 - 200 MeV) both in QCD and in QCD+. We can only be sure that n(fi) is a 
monotonically increasing function, which follows from the requirement of thermodynamic 
stability. 

The behavior of n(n) again becomes calculable in the region of very large fi 3> Aqcd- In 
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Figure 1: Schematic dependence of the baryon charge density on the chemical potential at T = 
(a) in QCD (jaq « m N - 16 MeV) and (b) in QCD+ (fj, « m N - 8 MeV). 



that case, the Pauli exclusion principle forces the quarks to occupy ever higher momentum 
states, and, due to asymptotic freedom, the interaction of quarks near the Fermi surface 
is (logarithmically) weak. The baryon charge density is proportional to the volume of a 
Fermi sphere of radius /i/3, n(/i) ~ iV/(/i/3) 3 /(37r 2 ). At low temperatures, only quarks near 
the Fermi surface contribute to the Debye screening of the gauge fields. The square of the 
screening mass, m 2 D , is proportional to the area of the Fermi surface: m 2 D ~ g 2 ^ 2 - This 



means that color interactions are screened on lengths 0(l/g/j) = 0(^\n(n/ Aqcd) / /•*) • This 
motivates the conclusion that nonperturbative phenomena such as chiral symmetry breaking 
should be absent at sufficiently large fi. Therefore, in QCD with massless quarks one should 
expect at least one other phase transition, at a value of /i which we define as /zi — a transition 
characterized by the restoration of the chiral symmetry. 

What is the value of /zi in QCD, and is it finite? Very little reliable information about 
the phase transition at /ii is available. However, several different approaches agree on the 
conclusion that the value of /ii is finite and that /xi — /i is on the order of the typical QCD 
scale Aqcd ~ 200 MeV ~ lfin -1 . For example, equating the quark pressure minus the MIT 
bag constant to the pressure of nuclear matter yields such an estimate (see, e.g., p0|). Here, 
we should also point out another interesting distinction between QCD and QCD+: the effect 
of the strange quark in QCD+ is to decrease the value of /ii compared to that of QCD. It 
has even been conjectured that this effect might be sufficient to drive /ii below /i , which 
would make normal nuclear matter metastable [IB], Another model which predicts the 
phase transition at finite fii is the Nambu-Jona-Lasinio model, which focuses on the degrees 
of freedom associated with the spontaneous chiral symmetry breaking and leads to a similar 
estimate for Hi 0. 

What is the order of this phase transition? The MIT bag model predicts that it is a first- 
order transition since the density, n, of the baryon charge is discontinuous. Unfortunately, 
analysis of the Nambu-Jona-Lasinio model shows that the order of the transition depends on 
the values of parameters, most notably, on the value of the cutoff. A larger cutoff leads to a 
second-order transition, a smaller to a first-order transition ||. A random matrix model at 
T = predicts a first-order phase transition M. In this paper, we shall extend the random 
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matrix model to permit consideration of the entire T\i plane. Before doing this, we shall use 
more general methods to analyze features of the phase diagram of QCD at finite density and 
temperature in the next section. 

An additional, qualitative argument for the first-order nature of the chiral phase tran- 
sition at Hi can be also drawn from a certain analogy of QCD to a metamagnet such as a 
crystal of ferrous chloride FeCl 2 . At temperatures below the Neel temperature, T N , and at 
zero magnetic field, H, such a crystal is antiferromagnetically ordered (i.e., the staggered 
magnetization, <p 8t , has a nonzero expectation value: (<p st ) 7^ 0. Analogously, (ipip) 7^ in 
QCD below T c . The magnetic field H is not an ordering field for the staggered magnetiza- 
tion because it couples to a different order parameter (i.e., normal magnetization, 0, with 
AS = —Hip) and induces nonzero (0). Similarly, the chemical potential induces nonzero 
('07o'i/'), and the term fit/jjoifj does not introduce explicit breaking of the chiral symmetry. 
At some critical value of H, ferrous chloride undergoes a first-order phase transition, and 
the staggered magnetization vanishes: (<p st ) = 0. One could naturally expect that in QCD a 
similar competition between the low temperature spontaneous ordering, {ipip) 7^ 0, and the 
ordering (^70^) 7^ induced by fi would result in a first-order phase transition. This anal- 
ogy can be continued into the Tfi plane or the TH plane in the case of the antiferro magnet. 
The antiferromagnet has a well known tricritical point in this plane. Its analogue in QCD 
will be discussed in section |[ 

Following the arguments of the two preceding paragraphs, we base our subsequent anal- 
ysis of the phase diagram of QCD with two massless quarks on the following expectations: 
(i) /ii ~ /io + 0(200 MeV) and (ii) the transition is of first order. 



4 Finite T and \i 

We shall use two order parameters to analyze the phase diagram of QCD at nonzero T and 
fi: the chiral condensate (ipip) (per flavor) given by 

^•)v--~ do) 

and the density of the baryon charge n given by eqn. (|j). We have already used n to show 
that there is a singularity at fi = fi^ and T = 0. It was important for that argument that n 
is exactly zero for all fi < fiQ. At nonzero T, however, n is not strictly for any fi > 0. For 
example, for very small fi and T one finds a very dilute gas of light mesons, nucleons and 
antinucleons with 

K( r,,) K £(^) 3/ V- f. (ii) 

Nevertheless, we can use a continuity argument to deduce that the first-order phase 
transition at T = 0, fi = /io has to remain a first-order phase transition for sufficiently 
small T. Therefore, there must be a line emerging from the point T = 0, fi = fiQ. One can 
think of this transition as boiling the nuclear fluid. The slope of this line can be related to 
the discontinuities in the entropy density, As (or the latent heat per volume TAs), and in 
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the baryon density, An, across the phase transition line through the generalized Clapeyron- 
Clausius relation: 

dfx As v ' 

This relation follows from the condition that the pressure, temperature and chemical po- 
tential should be the same in the two phases on a phase coexistence curve and eqn. (|). In 
analogy with ordinary liquid-gas transitions, the gaseous phase has a lower particle density 
(whence An < 0) and lower entropy density^ (whence As < 0). Therefore, the slope dT/dfi 
must be negative. We further expect that the slope is infinite at T = since s(T = 0) = 0, 
and hence As(T = 0) = 0. As there is no symmetry-breaking order parameter which dis- 
tinguishes the two phases, there is no reason why these two phases cannot be connected 
analytically. As in a typical liquid-gas transition, it is natural to expect that the first-order 
phase transition line terminates at a critical point with the critical exponents of the three- 
dimensional Ising model.0 The temperature of this critical point can be estimated from the 
binding energy per nucleon in cold nuclear matter, T = C(10MeV). (See Fig. |.) Signa- 
tures of this point are seen in heavy ion collisions at moderate energies (i.e., ~ 1 GeV per 
nucleon), and the critical properties of this point have been studied through measurements 



of the yields of nuclear fragments |2(| |2T|. In particular, the reported critical exponents are 
in agreement with those of the three-dimensional Ising model |2l| . 

Additional phase transitions which might occur at T = would give rise to additional 
phase transition lines. One could expect two generic situations. If there is a breaking of a 
global symmetry (e.g., translational symmetry in the case of nuclear matter crystallization), 
the phase transition line must separate such a phase from the symmetric phase at higher 
temperature without any gaps in the line. Otherwise, the transition can terminate at a 
critical point. 

At very high T ^> Aqcdj we have a plasma of quarks and gluons with a logarithmically 
small effective coupling constant, g(T), and we can again calculate the density of the baryon 
charge n: 



n(T, jj) 



d 3 p 
(2^ 



cxp 



\P\ ~/V3 
T 



n -l 



- {M -> -M} • (13) 



We expect that the chiral condensate is zero at very high T since the effective coupling is 
weak because of asymptotic freedom. Therefore, a phase transition must separate the quark 
gluon phase from the low temperature phase. This transition has been studied extensively 
at /i = using a variety of methods. In particular, lattice calculations have established the 
value of T c as approximately 160 MeV Arguments based on universality suggest that 
this transition is of second order with critical exponents of the SU(2)^x SU(2)# ~ 0(4) 
universality class [22]. Lattice calculations seem to confirm this scenario ^3|.f] Here, we 



lr Tlie entropy per particle is greater in the gaseous phase, but the entropy per volume s is smaller because 
of much smaller particle density. 

2 The Ising nature of the universality class follows from the fact that the transition can be modeled by 
an Ising lattice gas. 

3 A sufficiently light third quark would drive the transition first-order |22j . However, lattice calculations 
also indicate that the strange quark is not sufficiently light for this to occur [E4|. 
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assume that this is the case and try to understand what happens to this transition when \i 
is not zero. 

For massless quarks, the low-temperature hadronic phase and the quark-gluon plasma 
phase can be distinguished by the expectation value of (ipip), since this is identically zero in 
the quark-gluon phase and nonzero in the hadronic phase with spontaneously broken chiral 
symmetry. Therefore, when quark masses are strictly zero, a phase transition must separate 
these two phases, i.e., these phases cannot be connected analytically in the Tfi plane at 
m — 0. Therefore, a line of phase transitions must begin from the point T = T c , fi = and 
continue into the T\x plane. 

As discussed above, chiral symmetry restoration at T = is most likely to proceed via a 
first-order phase transition. Therefore, the transition must remain first-order as we continue 
along a line into the T/z plane. The slope of this line can again be related to the discontinuity 
in the baryon charge and the entropy density (|r2|) . Since we expect that both density and 
entropy will be larger in the quark-gluon phase, the slope of this line, dT/dfi, should be 
negative. 

This first-order transition line cannot terminate because the order parameter, (ipip), is 
identically zero on the one side of the transition. The minimal possibility is that it merges 
with the second-order phase transition line coming from T = T c , fi = 0; the point where the 
two lines join is a tricritical point || [7|, ||. Such a point exists in many physical systems 
(e.g., in the FeCk> antiferromagnet), and universal behavior in the vicinity of this point has 
been studied extensively. In the next section, we review those properties of a tricritical point 
which follow from universality. 

5 Universal properties of the tricritical point 

By analogy with an ordinary (bi) critical point, where two distinct coexisting phases become 
identical, one can define the tricritical point as a point where three coexisting phases become 
identical simultaneously. A tricritical point marks an end-point of three-phase coexistence. 
In order to see this in QCD, it is necessary to consider another dimension in the space of 
parameters — the quark mass m. This parameter breaks chiral symmetry explicitly. In such 
a three-dimensional space of parameters, one can see that there are two surfaces (symmetric 
with respect to m — > — m reflection) of first-order phase transitions emanating from the 
first-order line at m — 0. On these surfaces or wings with m ^ 0, two phases coexist: a 
low density phase and a high density phase. There is no symmetry distinguishing these two 
phases since chiral symmetry is explicitly broken when m/fl. Therefore, the surfaces can 
have an edge which is a line of critical points. These lines, or wing lines, emanate from the 
tricritical point. The first-order phase transition line can now be recognized as a line where 
three phases coexist: the high T and density phase and two low density and T phases with 
opposite signs of m and, hence, also of (ipip)- This line is called, therefore, a triple line. 

The plane m = is a symmetry plane. Chiral symmetry is exact only in this plane, 
and it is only here that the low and the high temperature phases must be separated by a 
transition. One can also view this plane as a first-order phase transition surface, since (t^ip) 
has a discontinuity across it. Then, the second-order phase transition line together with the 
triple line provide a boundary for this surface. 
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Critical behavior near the tricritical point can now be inferred from universality. The 
upper critical dimension for this point is 3. Since critical fluctuations are effectively three- 
dimensional for the second-order phase transition at finite T, we conclude that behavior 
near this point is described by mean field exponents with only logarithmic corrections. The 
effective Landau-Ginsburg theory for the long- wavelength modes, ~ ("0?/>), near this point 
requires a 6 potential which has the form (in the symmetry plane m — 0) 

fieff = n (T, n) + a(T, /i)0 2 + b(T, u)0 4 + c(T, /i)0 6 (14) 

with c > 0. The 6 term is necessary in order to create three minima corresponding to 
the three coexisting phases. This explains why the critical dimensionality is 3, since for 
this dimension, the operator 6 becomes a marginal operator. When b > 0, the transition 
occurs when a = and is a second-order transition similar to that seen in a 4 theory. This 
corresponds to the second-order line. When b < the transition occurs at some positive 
value of a and is of first order. This is the triple line. When both a and b vanish, we have a 
tricritical point. 

In particular, the following exponents in the symmetry plane m = are readily found 
using mean field 6 theory (as noted above, renormalization group studies show that 



the actual singularities include additional, logarithmic corrections). The discontinuity in the 
order parameter (0) = (ipip) along the triple line as a function of the distance from the 
critical point /i 3 , T 3 (measured either as T 3 — T or /i — /x 3 ) behaves like 

A(^)~(^-M 3 ) 1/2 . (15) 
The discontinuity in the density, n = dQ e g/dfi, across the triple line behaves like 

A"-^-^) 1 . (16) 

The critical behavior along the second-order line is everywhere the same as at the point 
fi = 0, T = T c (which is an infrared attractive fixed point). Therefore, (ipip) vanishes on the 
second-order line with 0(4) exponents. At the tricritical point, however, the exponent with 
which (ipip) vanishes is given by the Landau-Ginzburg theory as 

(^•0) ~ (T 3 — T) 1 / 4 . (17) 

When m/0, the potential f2 e fr(0) can also contain terms <fi and 3 which break — > — 
symmetry explicitly. (The term 5 can be absorbed by a shift of 0.) The potential (0) still 
has three minima, and a first-order phase transition can occur when two adjacent minima are 
equally deep. These transitions form a surface of first-order phase transitions — the wings. 
The two minima (and an intermediate maximum) can also fuse into a single minimum. This 
happens on the wing lines at the edge of the surface of first-order phase transitions. The 
critical behavior along the wing lines is given by the three-dimensional Ising exponents, as 
is usual at the endpoints of first-order liquid-gas type phase transitions not associated with 
restoration of a symmetry. In particular, the discontinuity in A(ifjtp) and An vanishes with 
exponent (3 ~ 0.31. These discontinuities are related to the slope of the wing surface at 
constant T through a relation similar to (p!2|): 

dfi A(i>ij)N f 



dm An 



(18) 
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There are many other universal properties in the vicinity of a tricritical point which can 
be derived from the above 6 Landau- Ginzburg effective potential. One can, for example, 
show that the m = second-order line, the wing lines, and the triple lines approach the triple 
point with the same tangential direction: The second-order line approaches from one side 
while the wing lines and the triple line approach from the opposite side. For more detailed 



description of the properties of tricritical points, see ref. [25 



6 A random matrix model at finite T and fi 

Random matrix models have proven to be a valuable tool for studying spontaneous chiral 
symmetry breaking in QCD. For example, it has been conjectured that the distribution of the 
eigenvalues of the Dirac operator near zero is universal [p6| . The universal expressions show 



a remarkable agreement with lattice Monte Carlo data |27| and are consistent with spectral 



sum rules from chiral perturbation theory |26| . Random matrix models have also been used 



to study chiral symmetry restoration phenomenon at finite temperature [28, 2£, 3C] as well 
as finite chemical potential [f|, [H], |33|] . 

Random matrix theory provides an effective description of those degrees of freedom in 
QCD which are responsible for the spontaneous breaking of chiral symmetry. In this respect, 
it is similar to Landau- Ginzburg effective theory. Random matrix theory is based on the 
observation that the spontaneous chiral symmetry breaking is related to the density of small 
eigenvalues of the Dirac operator (A <C Aq CD ). This relationship is expressed quantitatively 
by the Banks-Casher formula |34j], {ipip) = 7rp ev (0). Here, p ev (0) is the density of small 
(but non-zero) eigenvalues (per unit A and per unit four-volume, V4) of the Euclidean Dirac 
operator in the thermodynamic limit V4 —>■ 00. The dynamics of these eigenvalues can 
be described using a random matrix (of infinite size) in place of the Dirac operator. This 



approximation can be shown to give exact results in the mesoscopic limit |26|, |35 |. 

When the chemical potential is non-zero, the Dirac operator is not hermitian, and its 
determinant is no longer real. As a result, the density of its eigenvalues can be defined 
straightforwardly only in quenched QCD, i.e., when the contribution of the (complex) fermion 
determinant to the measure is approximated by unity. Fortunately, a more general relation 
exists between the chiral condensate and the linear density of zeroes of the partition function: 

(W) = vrp(O) . (19) 



As observed by Yang and Lee |3£|, non-analytic behavior in a thermodynamic quantity, 



including in the present case the discontinuity in the value of an order parameter in the 
thermodynamic limit, is caused by the coalescence of zeros of the partition function to form 
a boundary crossing the relevant parameter axis. In the case of QCD, the signature of 
spontaneous chiral symmetry breaking is the discontinuity in (iptp) as m is varied along the 
real axis and crosses m = 0. This discontinuity is equal to 27rp(0) where p(0) is the density 
of the zeroes on the imaginary m axis near m = in the thermodynamic limit. One can also 
show that p — > pcv in the quenched limit Nf — > but only when p = 0. For nonzero p, the 
density of the eigenvalues, p ev , is fundamentally different from the Nf — > limit of p 0). 

It can be shown that certain properties of the small eigenvalues of the Dirac operator are 
universal and are identical to those of a random matrix model with appropriate symmetries 
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37| , |38| , [35] , |27| , f40f . The virtue of the random matrix model is that it is solvable, 
i.e., one can calculate the distribution of the eigenvalues and of the Yang-Lee zeroes. The 
partition function of QCD at finite temperature and chemical potential in random matrix 
approximation is given by 

Z RM = J VXexp TrXxA det N f(D + m) , (20) 
where D is the 2N x 2N matrix approximating the Dirac operator p + (h/N^q. 



The random matrix X has dimension N x N . The total dimension of D is 2N. This 
is the number of small eigenvalues, which is proportional to V4. In QCD we expect iV 
to be approximately equal to the typical number of instantons (or anti-instantons) in V4; 
therefore, N/V4 ~ ni nst ~ 0.5 fm -4 |4l[]. The matrix C is deterministic and describes the 
effects of temperature and chemical potential. In the simplest (and original) T 7^ 0, /J = 



model [28], the choice C = ttT describes the effect of the smallest Matsubara frequency. As 



noted in ref. pOfl , it is possible to simulate the effects of the eigenvalue correlations induced 
by the pairing of instantons and anti-instantons into molecules by choosing a more general 
form for the diagonal matrix C with elements, Ck, which are (increasing) functions of T. 
In the T = 0, /i 7^ model of ref. ||], C — n/{iN c ) describes the effect of the chemical 
potential. In this paper, we consider the more general case T 7^ 0, /x 7^ 0. Although we 
do not know the detailed dependence of the elements of C on T and /i, we understand that 
T primarily affects the real (i.e., hermitian) part of C and \i affects the imaginary (i.e., 
antihermitian) part. We shall adopt the following approximate form for this dependence, 
Ck = cnrT + bfi/(iN c ) for one half of eigenvalues and Ck = — airT + bn/ (iN c ) for the other half 
with a and b dimensionless parameters.^ This form accounts for the fact that there are two 
smallest Matsubara frequencies^], +ttT and — irT . Such a linear ansatz for C is certainly very 
naive, but in this paper we decided not to try and refine it. This form reflects sufficiently 
well our understanding of the properties of C. 
The chiral condensate is calculated as: 

1 7 /\ 1 <91nZ RM , x 

m = n^-^t- (22) 

Current algebra fixes the value of (ipip)o ~ 2fm -3 at T — fj, — m — 0. The only dimensionful 
parameter remaining in the partition function, Zrm, is the variance of the random matrix, 
a. Thus, 

('ipip)o = const W const — ^ . (23) 

V4CX a 



4 These parameters are intended to reflect the degree of overlap and correlation between instantons and 
anti-instantons. One can therefore anticipate that a and b are smaller than 1. We shall estimate the values 
of a and b below. 

5 Alternatively, this form preserves the relation (70-D) = at fj, = 0. 
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The dimensionless constant will be found below and is equal to 2. This fixes the value of 
cr « 0.5 fm _1 rs 100 MeV. It is convenient to use a as a unit of mass in the model and also 
absorb the coefficients na and b/N c into T and [l. In other words, we measure m in units of 
cr, T in units of a/(ira) and // in units of aN c /b. 

The N —>■ oo (i.e., thermodynamic) limit of the partition function, eqn. (pip, can be 
found in the now-standard way |2£| |29], ^0| . Performing the Gaussian integration over X 
and introducing auxiliary Nf x Afy matrices 0, one can rewrite the partition function in the 
form 



jv*~*-™mw ( IXTt 7:1 ) ^ m ( r.™ *r + ™ ) 

£>0exp[-iVft(0)], (24) 



where 



ft(0) = Tr[00 f - - ln{[(0 + m)(0 f + m) - (/i + iT) 2 ] • [(0 + m)(0 f + m) - (// - iT) 2 ]}]. (25) 
2 

The integration in eqn. (|24]) is performed over 2 x Nf x iVy variables which are the real and 
imaginary parts of the elements of the complex matrix 0. In the limit N —>■ oo this integral 
is determined by a saddle point of the integrand or, alternatively, the minimum of Q(<f)): 

lim — In Z RM = - min fi(0) . (26) 

The function Q(4>) is an effective potential for the degrees of freedom describing the dynamics 
of the chiral phase transition. One can see that the value of <fi at the minimum, i.e., the 



equilibrium value (0) , gives us the value of the chiral condensate, eqn. (|22|) : 



<^> = -±--2ReTr(0>, (27) 
ivy 1/4 a 

(cf. eqn. (|23|) ). 

For real m, it is reasonable to expect that the minimum occurs when is a real matrix 
proportional to a unit matrix. With this assumption, we need to find only one real parameter, 
0, which minimizes the potential: 



Q = N 



\ In {[(0 + m) 2 - (/x + iT) 2 } • [(0 + m) 2 - {jjl - iTf]} 



(28) 



This is not a simple 6 potential, but it has remarkably similar properties. In particular, the 
condition dVL/dcj) = gives a fifth-order polynomial equation in 0. Let us first consider the 
symmetry plane m = 0. Then the equation dQ/d<f) = always has one trivial root, = 0. 
The remaining four roots are the solutions of a quartic equation, which has the form of a 
quadratic equation in 2 : 

a4 O ( ..2 rp2 1 1\ i2 1 I .2 1 rp2\2 1 .2 



_ 2 ( ^ - T 2 + - ) 2 + (/i 2 + T 2 ) 2 + /i 2 - T 2 = 0. (29) 
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Above the second-order line, i.e., in the high temperature phase, eqn. ( p9| ) does not have real 
roots (since 2 < for each pair of roots). This corresponds to the fact that the potential 
Q has only one minimum at = (i.e., the trivial root). On the second-order line, a pair 
of roots of eqn. (|29|) becomes zero, i.e., the potential is Q ~ 4 near the origin. This means 
that on the second-order line, 

(/i 2 + T 2 ) 2 + /x 2 -T 2 = . (30) 

The second-order line ends when the remaining pair of roots also becomes zero, i.e., the 
potential becomes Q ~ 6 near the origin. This happens when 

/i 2 -T 2 + ^ = (31) 

on the second-order line. The condition, eqn. (|31"|) together with eqn. (|30D determines the 
location of the tricritical point in the T\x plane: 

T 3 = ^vV2 + 1 w 0.776 and /j 3 = ^y/ y/2 - 1 « 0.322 . (32) 

The equation for the triple line is obtained from the requirement that the depth of the 
minima in Q at given by the pair of solutions of eqn. (p9|) (farthest from the origin) should 
coincide with the depth at the origin = 0. The equation for the triple line is therefore 

f - T 2 + \ + iy/l-ie^T*- \ In ^ + V1~16^ + ln(//2 + T2) = Q (33) 

In particular, when T = 0, we obtain the elementary equation /i 2 + 1 + In /i 2 = whose 
solution is fi ~ 0.528 ||. This is the value of H\. Setting \i = in the equation for the 
second-order line, eqn. (|30|) , we find T c = 1 |28| . 

We recall that the units of T and /x depend on the unknown dimensionless parameters a 
and b. However, these unknown factors cancel from the ratios 

T 

« 0.78; — w 0.61. (34) 

Taking T c = 160 MeV and /i X = 1200 MeV, we find that T 3 w 120 MeV and /i 3 w 700 MeV.Q 
Note that the second-order line, eqn. (|30D, marks the location of the points on the phase 
diagram where the symmetric minimum = disappears, i.e., turns into a maximum. 
Continuing this line below point T 3 , we obtain the location of spinodal points. In the region 
between this line of spinodal points and the first-order phase transition line, the chirally 
symmetric phase = can exist as a metastable state. Such a state can be reached 
by supercooling, and it is unstable towards the nucleation of bubbles of the broken phase 
7^ 0. A similar line with equation 4T/x = 1, the superheating line, together with the 
supercooling line, eqn. fl30D , bound the region around the first-order phase transition line 
where the potential Q(4>) has 3 minima. All these lines meet at the tricritical point. 

6 With these values for T c and fii, the parameters a and b have the values a w 0.2 and b w 0.13. 



14 



Away from the symmetry plane m = 0, the expressions for the wing surfaces and the 
wing lines become rather lengthy and will be presented elsewhere. The principle, however, 
remains simple. The minima of the potential, Q, satisfy the equation dQ/d<f) = and are 
given by (three out of five) roots of a fifth-order polynomial. On the wing surface, the depth 
i n a P a h °f adjacent minima is the same. On the wing line, the two adjacent minima 
fuse into one. In terms of the roots of the polynomial, three roots coincide (two minima and 
one maximum). In other words, the potential is Q ~ (<f) — (<p)) A on the wing line and near 
the minimum. 



T 




Figure 2: Phase diagram of QCD with two light flavors of mass m as calculated from the random 
matrix model. The almost parallel curves on the wing surface are cross sections of this surface with 
m =const planes. The units of m are a ~ 100 MeV, of T are T c ~ 160 MeV, of fi are //i/0.53 ~ 2300 
MeV, with the choices of T c and \i\ from the text. 

The resulting phase diagram is plotted in Fig. |2|. One can see, as expected from mean field 
theory near the tricritical point, that the wing lines together with the triple line approach 
the tricritical point with the same slope as the second-order line but from the other side. 
The critical exponents near the tricritical point as given by the random matrix model can 
be also seen to coincide, as expected, with the mean field exponents of eqns. (|l|), (|16D, and 

From the random matrix model, we also learn how the zeros of the partition function 
in the complex m plane evolve with changes in temperature and chemical potential. A few 
typical cases are illustrated in Fig. |3|. At zero T and fi, the zeros form a cut (in the iV — ► oo 
limit) along imaginary axis. Raising the temperature pushes the zeros away from the origin 
along the imaginary axis until the density at the origin vanishes (continuously), the cut 
breaks in two, as in Fig. |3]a, and chiral symmetry is restored (cf. eqn. fll9|)). The chemical 
potential pushes the zeros away from the origin in the direction of the real axis until the cut 
splits in two, as illustrated in Fig. |b. Note that the density p(0) is finite just before the 
split. Therefore, the transition is of first order. Near the tricritical point, the split in the 
direction of the real axis (due to the chemical potential) occurs at the same time that the 
density p(0) vanishes (due to the effects of the temperature). This is illustrated in Fig. ^]c. 
The resulting dependence of (ipip) on T and /x is shown in Fig. |j. 

A comment should be added regarding the calculation of the baryon density, n, in the 
random matrix model. The value of (l/V^dln Z^_jd\x does not represent the complete 



15 



-1 



-2 — 



-3 



N=40 



T=1.000~ 
jU-0.000 



2 — 



-1 



-2 



-3 



N=40 



T=0.000 
/i=0.555 - 



- 3 - M= 



N=40 



— 2 — 



1 - 



— — 



2 — 



3 — 



T=0.78CT 
jU-0.390 



(a) (b) (c) 

Figure 3: Zeros of the partition function of a finite size N random matrix model ( |20| ) in the 
complex m plane calculated numerically at different values of T and \x. The calculation is done for 
Nf = 1, but the N — * oo limit is Nf independent. The density of points is proportional to the 
strength of the cut (discontinuity in {tpijj}) in the N — » oo limit. 



baryon density. The reason is that In Z-rm contains only contributions from the soft modes 
of the condensate, <fi ~ ipip. Further dependence on \x is contained in the contributions to 
In Z from other degrees of freedom. In the effective Landau- Ginzburg theory, such additional 
contributions are embedded in the term fl (T,fj,), eqn. (0). However, the singular behavior 
of the system is exclusively due to the soft modes of the condensate, i.e., terms involving 
in Q eS . Thus, it is legitimate to calculate singular properties of n, such as An, near a 
(tri)critical point using (1/ V^jdln Z RM /dfi. At T = 0, for example, we find: 

= NN f (Jj- + w 5NNf, (35) 

in units of which is 6/A^ c S. With our previous choice of jii = 1200 MeV (i.e., b w 0.13), 
we find that An w 0.4 fm^ 3 « 2.5 n , which seems a reasonable estimate. 



AnV, = N 




<t>=o 




7 Discussion and summary 

In this paper, we have presented an analysis (qualitative and, in some cases, quantitative) 
of the salient features of the phase diagram of QCD with two light or massless quark flavors 
at finite temperature and baryon chemical potential. The most important features of this 
phase diagram are summarized in Fig. [|. The phase diagram can certainly have a much richer 
structure. The phase transitions shown there are distinguished by the fact that a good order 
parameter can be associated with each of them. Here the term "good order parameter" 
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Figure 4: The chiral condensate (ipi/)) (in units of (ipi/j)o ~ 2fm 3 ) as a function of T and /i in the 
random matrix model. The units of T and fj, are as in Fig. ||. 

implies the existence of some quantity whose expectation value is identically zero in some 
finite region of parameter space or in one phase and is some function of parameters in the 
other phase. Two such phases must be separated by a nonanalytic boundary, i.e. a phase 
transition. What is crucial here is the identical vanishing of an order parameter or its strict 
independence of the parameters of the theory. Usually, this is ensured by the existence of 
some symmetry with respect to which this order parameter transforms nontrivially. 

For the chiral phase transition, a good order parameter is the value of (ipip), which 
spontaneously breaks the global SU(2)^x SU(2) fi chiral symmetry to SU{2) V . Hence, the 
phase with (ipip) = and the phase with (ipip) 7^ cannot be connected without crossing a 
phase transition line in the T\x plane. 

The transition from n = to n 7^ along the T = line provides another example of 
a phase transition associated with a good order parameter. The phases n = and n / 
cannot be analytically connected; i.e. one must pass through a nonanalytic boundary (phase 
transition) when passing from one to the other. Since this is a first-order phase transition, 
continuity requires that there is also a first-order transition line for some T ^ 0. This line 
can, however, terminate since n is no longer a good order parameter when T^O. 

The existence of a good order parameter is a sufficient but not a necessary condition for 
a phase transition. Other phase transition lines associated with more subtle phenomena are 
also possible. One interesting example, which attracted attention recently, is the transition 
associated with color superconductivity. The existence of a color superconducting phase was 
first argued for by Bailin and Love on the grounds that one gluon exchange is attractive in 



the color antitriplet quark-quark channel |14| . This means that the Fermi surface at very 



high fi becomes unstable and forms a gap. This phenomenon was recently reanalyzed using 
other methods at moderate values of fi with the conclusion that the effect is enhanced by 
instanton-induced interactions [1TB, [15]. This means that another finite-T transition, which 



stretches all the way to /i = 00, may be present on the phase diagram, Fig. [5| 

This transition, unlike the chiral phase transition and the nuclear matter liquid-gas tran- 
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M-o M-i 00 

Figure 5: A schematic phase diagram of QCD with 2 massless quark flavors. Other phase transition 
lines are possible, for example, in the low temperature region to the right of fiQ. Another example 
is a transition associated with color superconductivity plotted as a dashed line. Thicker lines are 
first-order phase transitions. The T c — T3 line is a second-order phase transition. The tricritical 
point is at T3, ^3 and the critical point of the nuclear matter liquid-gas transition is at To. 

sition (at T = 0), does not seem to have a good order parameter associated with it. In 
particular, the diquark condensate, (ipip), is not gauge invariant. The dynamical mechanism 
responsible for the binding of diquark pairs is certainly operative at low temperatures, but 
the absence of a good order parameter does not allow us to assert that the temperature 
induced transition associated with the breaking of diquark pairs must always (i.e., at all fi) 
proceed through a thermodynamic singularity rather than a smooth analytic crossover. 

For some purposes, it is more natural to study the phase diagram in the space of density 
and temperature. The phase diagram of Fig. [| can be converted into such a space and is 
shown in Fig. ||. A typical feature is that the first-order phase transition line from the T\x 
plane now appears as a region of phase coexistence in Fig. || In equilibrium, the values of 
density and temperature inside this region can be achieved only through an inhomogeneous 
mixture of two phases with different densities but the same T and /i. These densities are 
indicated by the ends of the horizontal lines drawn in the phase coexistence region. 

The most interesting feature of the phase diagram of Figs. |5| and ^] is the presence of 
a tricritical point. Because of the fact that the critical dimensionality for such a point is 
equal to 3, critical behavior near this point is given by mean field theory plus logarithmic 
corrections. In particular, a simple random matrix model predicts the correct algebraic 
critical exponents. 

The tricritical point lies in the region expected to be probed by heavy ion collision 
experiments. It would be interesting to find an experimental signature for such a point. 
Since quark masses are not precisely zero, we should consider a slice of a three-dimensional 
phase diagram, Fig. ^|, with m 7^ 0. A qualitative difference between the phase diagram 
for m 7^ and that for m = is the absence of the second-order phase transition line 
associated with the restoration of chiral symmetry. This symmetry is explicitly broken for 
m^O. However, continuity from m = ensures that the first-order finite density transition 
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Figure 6: The phase diagram of Fig. || shown as a function of density and temperature. The 
dashed line of the color superconductivity transition is not drawn. Horizontal lines connect points 
corresponding to densities of phases on two sides of the first-order line (i.e., the coexistence curve) 
of Fig. [|. The points n$ on the T = line are the same as on Fig. [l]a. 



is still present at m ^ 0. This transition line is terminated by an ordinary critical point. 
Criticality at this point is not associated with chiral symmetry restoration, and excitations 
with the quantum numbers of pions do not become massless there. Criticality at this point is 
associated with the fact that a correlation length in the channel with the quantum numbers 
of the sigma meson becomes infinite. (Hence, it is plausible to infer that this point has the 
critical behavior of the three-dimensional Ising model.) Possible experimental signatures of 
this phenomenon are under investigation. 
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Note added 

After this work was completed, a paper [|2| appeared which addresses similar questions in 
the context of a Nambu-Jona-Lasinio model for color superconductivity. The results of [i2| 
agree with and complement our findings. 
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